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Magnetic properties of the spin-3/2 Heisenberg system Cr2[BP30i2] are investigated by magnetic 
susceptibility x(T) measurements, electron spin resonance, neutron diffraction, and density func- 
tional theory (DFT) calculations, as well as classical and quantum Monte Carlo (MC) simulations. 
The broad maximum of x(T) at 85 K and the antiferromagnetic Weiss temperature of 139 K indi- 
cate low-dimensional magnetic behavior. Below Tn = 28K, Cr2[BP30i2] is antiferromagnetically 
ordered with the k = propagation vector and an ordered moment of 2.5 /j,b/Ct. DFT calculations, 
including DFT+t/ and hybrid functionals, yield a microscopic model of spin chains with alternat- 
ing nearest-neighbor couplings J\ and J[. The chains are coupled by two inequivalent interchain 
exchanges of similar strength (~1-2K), but different sign (antiferromagnetic and ferromagnetic). 
The resulting spin lattice is quasi-one-dimensional and not frustrated. Quantum MC simulations 
show excellent agreement with the experimental data for the parameters Ji ~ 50 K and J[/ Ji~0.5. 
Therefore, Cr2[BP30i2] is close to the gapless critical point (J(/Ji = 0.41) of the spin-3/2 bond- 
alternating Heisenberg chain. The applicability limits of the classical approximation are addressed 
by quantum and classical MC simulations. Implications for a wide range of low-dimensional S — 3/2 
materials are discussed. 

PACS numbers: 75.50.Ee, 75.30. Et, 75.10.Hk, 75.10.Jm, 75.10.Pq, 61.66.Fn 



I. INTRODUCTION 

Magnetic properties of transition-metal compounds are 
generally described with the Heisenberg spin Hamilto- 
nian that may by augmented by additional terms re- 
sponsible for the anisotropy. For example, the following 
Hamiltonian 

i/ = ^JS 4 -S,-^A(S l -z) 2 (1) 

(ij) ' 

accounts for a variety of experimental situations. Here, 
the first term describes the isotropic coupling J between 
spins on sites i and j, the second term is the single- ion 
anisotropy Di, and z is a unitary vector along z, whereas 
Si and Sj are quantum-mechanical spin operators. Al- 
though a handful of solvable Heisenberg models exist, 1 ' 2 
the complex algebra of the spin operators generally im- 
pedes analytical solutions. 

A widely used approximation of a Heisenberg model is 
its classical treatment: the original quantum mechanical 
spin operators are replaced by real-space vectors <Sj. 
This transformation leads to enormous alleviation of the 
computational effort, even if the topology of exchange 
couplings is very intricate. 

A fundamental limitation of classical models is their 
inability to account for quantum-mechanical singlets and 
the ensuing underestimation of the ground state (GS) 
energy. Indeed, for an antiferromagnetic (AF) exchange 
coupling J, the energy of the quantum-mechanical singlet 



state —JS(S + 1) is always lower than the classical en- 
ergy — JS 2 of the antiparallel arrangement. At the same 
time, the relative energy gain of a singlet scales as S^ 1 , 
thus being maximal for S—l/2 and infinitesimal in the 
S — > oo limit. 

Presently, there is an empirical evidence that for S> 2, 
classical models capture the essential physics and cor- 
rectly reproduce the experimental magnetic behavior. 3-5 
In contrast, classical models often fail to predict the 
correct GS and magnetic excitation spectrum for the 
extreme "quantum" case of S = 1/2. 6,7 Further on, 
many 5=1 systems, e.g., one-dimensional (ID) Haldane 
chains, 8 can not be treated classically. Therefore, it is 
crucial to establish the applicability limits for the clas- 
sical approximation, in order to distinguish between the 
"classical" cases amenable to a simplified model treat- 
ment, and the "quantum" cases that require the complete 
quantum-mechanical solution of the spin Hamiltonian. 

A feasible way towards better understanding of these 
limits are real material studies, that allow for a direct 
comparison between theory and experiment. Previous 
work on the quantum-classical crossover rendered the 
quasi-lD S = 3/2 magnets as the relevant playground. 
Despite the relatively large spin of 3/2, the pronounced 
one-dimensionality may lead to sizable quantum fluctu- 
ations and, thus, to deviations from the classical behav- 
ior. For example, inelastic neutron scattering (INS) stud- 
ies of Cs[VA 3 ] (X = G1, Br) evidence that at high tem- 
peratures these materials are classical, while lowering 
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the temperature results in a crossover to the quantum 
behavior. 9 The INS experiments on the quasi-lD mag- 
net AgCr[P 2 S 6 ] yielded a sizable discrepancy between 
the observed spin-wave velocity and its classical value, 
thus indicating strong quantum effects that are present 
in S = 3/2 chains. 10 

Here, we present a joint experimental and theoreti- 
cal study of the quasi-lD S = 3/2 system Cr 2 [BP 3 0i2]. 
Its crystal structure features magnetic Cr 3 , -1- Og blocks 
embedded into a complex borophosphate framework 
(Fig. 1). Although this type of crystal structure could 
lead to the simple magnetism of isolated spin dimers, 
neutron diffraction and magnetic susceptibility measure- 
ments evidence long-range AF ordering at Tn = 28 K 
that is indicative of sizable interdimer couplings. Us- 
ing extensive density-functional theory (DFT) calcula- 
tions, we evaluate a microscopic magnetic model for this 
compound and establish the spin lattice of weakly cou- 
pled bond-alternating Heisenberg chains. This theoreti- 
cal model is in good agreement with the neutron scatter- 
ing and magnetic susceptibility data. 

As the microscopic magnetic model of Cr2[BPaOi2] 
lacks frustration, its properties can be simulated with the 
computationally efficient quantum Monte Carlo (QMC) 
techniques. At the same time, the classical model on the 
same spin lattice can be treated using classical Monte 
Carlo (MC) algorithms. By comparing the QMC and 
classical MC results, we evaluate the relative impact of 
quantum as well as thermal fluctuations on the spin cor- 
relations, and in this way address the crossover between 
the quantum and classical behavior. 

This paper is organized as follows. The experimen- 
tal as well as numerical methods are described in Sec. II. 
The experimental part in Sec. Ill comprises the results of 
neutron diffraction, electron spin resonance (ESR), and 
magnetic susceptibility measurements. The DFT-based 
evaluation of the microscopic magnetic model (Sec. IV) 
is followed by the refinement of model parameters by 
means of QMC simulations and subsequent fitting to the 
experiment (Sec. V). The differences between the quan- 
tum model and its classical counterpart are discussed in 
Sec. VI. We summarize our results and give a short out- 
look in Sec. VII. 



II. METHODS 

Cr 2 [BP 3 0i2] was prepared by the Pechini-type 
method. 11 2.172 g chromium acetate (Alfa Aesar, 
23.37(5) wt% chromium content) and 5.7945 g citric acid 
(Alfa Aesar 99+%) were dissolved in 50 ml water. 5 ml 
of glycerol (Sigma 99%) and a well-ground mixture of 
0.3197 g boric acid (Sigma) and 1.7613 g ammonium dihy- 
drogen phosphate (Merck) were added to the dark-green 
solution. After slow evaporation of water, a transparent 
resin was formed which was dried at 200 °C for 2h. The 
obtained product was crushed and transferred to a corun- 
dum crucible. The first heating at 850 °C for 12 h yielded 



a grey-green product. Several annealing steps with grind- 
ing in between followed. The final product was obtained 
after heat treatment at 1000 °C for 48 h. 

The powder X-ray diffraction pattern (Huber image 
plate Guinier camera G670, Ge monochromator, Cu 
K a i radiation, A = 1.5406 A, powdered sample fixed with 
Vaseline between two Mylar foils each 6 /jm thick) indi- 
cated the formation of single-phase Cr2[BP30i2]. 12 

Neutron powder diffraction data were collected 13 using 
the BT-1 high-resolution powder diffractometer at the 
NIST Center for Neutron Research. A monochromatic 
neutron beam with the wavelength of 1.5403 Awas pro- 
duced by a Cu (311) monochromator. Collimators with 
horizontal divergences of 15', 20', and 7' full width at 
half-maximum were used before and after the monochro- 
mator, and after the sample, respectively. The inten- 
sities were measured in steps of 0.05° in the 28 range 
3-168°. The data were collected at 4K, 35 K, and 300 K. 
Additionally, the magnetic scattering was studied with 
the triple-axis spectrometer BT-7 using the wavelength 
of 2.359 A. The intensity of the strongest magnetic reflec- 
tion was monitored with the step of 0.5 K in the 5 — 35 K 
temperature range. The structural analysis was per- 
formed using the programs GSAS. 14 The magnetic struc- 
ture was refined with Fullprof . 

The magnetic susceptibility was measured using a com- 
mercial Quantum Design MPMS SQUID in the tempera- 
ture range 2-380 K in magnetic fields up to 5 T. The elec- 
tron spin resonance (ESR) measurement was performed 
at room temperature with a standard continuous-wave 
spectrometer at X-band frequencies (y ~ 9.5 GHz) by 
using a cylindrical resonator in TE012 mode. 

DFT calculations have been performed using the 
full-potential local-orbital code fplo9. 00-33 (Ref. 16) 
and the pseudopotential projector-augmented-wave code 
vasp-5.2 (Ref. 17). For the scalar relativistic 
calculations, we used the local density approxima- 
tion (LDA) 18 and generalized gradient approxima- 
tion (GGA) 19 exchange-correlation potentials. Spin- 
unpolarized calculations were performed on a 14x14x14 
mesh of fc-points. For the spin-polarized calculations 
(DFT, DFT+J7), we doubled the cell along c and used a 
4x4x2 fc-mesh. Hybrid- functional calculations were per- 
formed in VASP using the HSE06 functional 20 on a 2x2x2 
fc-mesh. 21 The convergence with respect to the /c-meshes 
has been accurately checked. All calculations have been 
performed based on the crystal structure determined at 
4 K as given in Table I. 

QMC simulations were performed using the code 
looper from the software package ALPS version 1.3. 22 
The magnetic susceptibility was simulated on 8x8x64 
finite lattices of 5 = 3/2 spins in the temperature range 
0.25-8 Ji, corresponding to 12.5-400 K (see Sec. V) using 
30 000 sweeps for thermalization and 300 000 sweeps after 
thermalization. The statistical errors (<0.5%) are below 
the experimental accuracy. For simulations of the spin 
stiffness and the static structure factor, we used finite 
lattices up to 13 824 and 2 048 spins, respectively. 
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For the classical MC simulations, we used the SPINMC 
code, 22 with 200 000 and 2 000 000 sweeps for and after 
thermalization, respectively. Chains of N = 800 spins 
were evaluated. Exact diagonalization of the S z =0 sec- 
tor was performed for N = 14 sites alternating S = 3/2 
chain using SPARSEDIAG from the ALPS package. 22 



III. EXPERIMENTAL RESULTS 
A. Crystal structure 

The crystal structure of Cr2[BP30i2] has been ini- 
tially refined from X-ray powder data in the space 
group P3. 12 However, the crystal structures of related 
M(III) 2 [BP 3 0i2] borophosphates were determined (M = 
In) or re-determined (M = Fe) from single crystal data 
in the space group P6 3 /m. 23 This apparent mismatch led 
us to reconsider the crystal structure of Cr2[BP30i2]- A 
thorough analysis of both X-ray and neutron powder pat- 
terns identified the reflection condition 00Z with I = In, 
which is characteristic of the 63 screw axis. No other re- 
flection conditions could be observed, so that the list of 
possible space groups is restricted to P63, PG^/m, and 
P6 3 22. While the refinement in P6322 was unsuccess- 
ful, the space group P6 3 /m resulted in low residuals and 
a fully ordered crystal structure. Therefore, this cen- 
trosymmetric space group was preferred over its non- 
centrosymmetric subgroup P6 3 . No significant struc- 
tural changes between 4K and room temperature were 
detected. Crystallographic data are listed in Table I. Se- 
lected interatomic distances at 4 K are given in Table II. 

The structure solution of Cr 2 [BP 3 0i2] in the space 
group P63/TO is consistent with earlier results for other 
7\/ 2 [BP 3 0i2] transition-metal borophosphates. 23 Their 
crystal structures are isotypic and feature dimers of face- 
sharing MOs octahedra (Fig. 1, top). The octahedra ex- 
hibit a sizable trigonal distortion leading to two groups of 
nonequivalent Cr-0 distances. The shorter Cr-Ol dis- 
tances of about 1.91 A (at 4K) take the terminal position 
of the C^Og dimer, whereas the longer Cr-03 distances 
of about 2.04A are bridging (Fig. 1, bottom right). The 
formation of longer Cr-0 distances in the bridging po- 
sition is due to the Cr-Cr contact remaining relatively 
long (about 2.8 A), in order to reduce the repulsion be- 
tween the positively charged Cr 3+ ions. PO4 tetrahedra 
link the Cr 2 Og dimers along [0001], whereas the BO3 tri- 
angles interconnect the PO4 tetrahedra and do not share 
oxygen atoms with the CrC>6 octahedra (Fig. 1, top). 



B. Magnetic structure 

At 4K, an additional magnetic scattering was ob- 
served. The only visible magnetic reflection is at 29 ~ 
17. 6° and matches the weak 101 reflection of the nuclear 
structure. As the Neel temperature of Cr2[BP30i 2 ] is 




FIG. 1. (Color online) Crystal structure of Cr 2 [BP 3 Oi2]. Top: 
view along [0001]. Bottom left panel: view perpendicular 
to [0001]. The pathways of the leading interchain couplings 
Jid and Ji C 2 are shown with dark gray (dark blue) and light 
gray (red) lines, respectively. Thick arrows denote the exper- 
imental magnetic structure. Bottom right panel: dimers are 
connected via three PO4 tetrahedra and form chains running 
along [0001]. The intra-dimer coupling Ji (double line) as 
well as the interdimer coupling J[ (single line) are shown. 



Tjv — 28 K (see below) and no structural changes be- 
low Tjy are expected, the subtraction of the 35 K data 
from the 4K data results in the purely magnetic scat- 
tering (Fig. 2, inset). However, no features other than 
the 101 magnetic reflection could be found. This solitary 
magnetic reflection is naturally indexed with the k = 
propagation vector which together with the 4/ position 
of Cr allows for several irreducible representations. The 
magnetic structure was successfully refined using the Ti 
representation that entails Cr spins directed along the c 
axis, as well as the antiparallel arrangement of the neigh- 
boring spins both in the ab plane and along the c direc- 
tion. The resulting magnetic structure is shown in Fig. 1, 
bottom left. 

To evaluate the ordered magnetic moment, we refined 
the subtracted pattern as the purely magnetic phase, and 
fixed the scale factor according to the refinement of the 
35 K pattern, which is free from magnetic scattering. The 
resulting magnetic moment at 4 K is 2.5(1) /is, slightly 
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TABLE I. Crystallographic data for Cr2[BP30i2] (space group PQ^/m) according to the neutron powder diffraction data at 4, 
35 and 300 K. 



T 




4K 






35 K 






300 K 




a (A) 




7.9444(2) 






7.9448(2) 






7.9524(2) 




C (A) 




7.3439(3) 






7.3448(3) 






7.3543(3) 




v (A 3 ) 




401.40(2) 






401.49(3) 






402.78(3) 




R P (%) 




5.5 






5.8 






6.4 




^wp (%) 




6.9 






7.6 






7.9 




-Rcxp (%) 




6.1 






7.5 






8.4 






x/a 


y/b 


z/c 


x/a 


y/b 


z/c 


x/a 


y/b 


z/c 


Cr 4/ 


1/3 


2/3 


0.0578(8) 


1/3 


2/3 


0.0614(9) 


1/3 


2/3 


0.0590(7) 


P 6h 


0.3657(4) 


0.3166(5) 


1/4 


0.3658(4) 


0.3168(5) 


1/4 


0.3653(4) 


0.3159(5) 


1/4 


B la 








1/4 








1/4 








1/4 


Ol 6h 


0.4023(4) 


0.5254(4) 


1/4 


0.4009(4) 


0.5249(5) 


1/4 


0.3997(4) 


0.5255(5) 


1/4 


02 6h 


0.1362(4) 


0.1931(4) 


1/4 


0.1351(5) 


0.1924(4) 


1/4 


0.1358(4) 


0.1926(4) 


1/4 


03 12i 


0.4368(3) 


0.2691(3) 


0.0774(2) 


0.4365(3) 


0.2697(4) 


0.0768(3) 


0.4382(3) 


0.2716(3) 


0.0786(3) 



TABLE II. Selected interatomic distances (in A) in the 
Cr2[BP30i2] crystal structure at 4K. 



Cr-OI 
Cr-03 
Cr-Cr 
B-02 



3 x 2.041(5) 
3 x 1.910(3) 

2.823(11) 
3 x 1.365(3) 



P-Ol 
P-02 
P-03 



1.535(4) 
1.580(4) 
2 x 1.511(2) 



lower than the spin-only value of 3 hb expected for the 
5 = 3/2 ion. Its temperature evolution can be tracked by 
the temperature dependence of the magnetic reflection 
101 (Fig. 2). These data are fitted with the empirical 
formula 



J(T) = 7, 



1 - 



T 



A' 







(2) 



where /bg refers to the nuclear scattering and back- 
ground above Tjv. The fit yields the Neel temperature 
T N = 28.2(7) K and the critical exponent f3 = 0.38(3). 
The estimated T/v — 28 K is in excellent agreement with 
the magnetic susceptibility measurement presented below 
(Fig. 3), whereas (3 falls into the range of values (0.36- 
0.39) proposed for the 3D Heisenberg model. 24 



C. Magnetic susceptibility and ESR 

The temperature dependence of the magnetic suscepti- 
bility x{T) (Fig. 3, top) reveals a typical low-dimensional 
behavior with a broad maximum around T — 85 K. The 
high-temperature part (T > 200 K) of the curve obeys 
the Curie-Weiss law (Fig. 3, bottom left) with the 
Curie constant C — 1.987 emu K (mol Cr) -1 leading to 
the effective magnetic moment ^ e ff = 3.987 fiB per Cr, 
which is slightly larger than the spin-only contribution 
fj, cS — 3.88 /is (assuming the orbital moment is corn- 
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FIG. 2. (Color online) Cr2[BPsOi2]: temperature evolution of 
the magnetic 101 reflection (circles) and the fit with Eq. (2). 
The residual intensity above Tjv is due to the nuclear scat- 
tering and background. Inset: refinement of the subtracted 
(Iak — ^35k) pattern. Note that the temperature scan (main 
figure) and the angular scan (inset) are done on different in- 
struments, hence the respective intensities should not be com- 
pared. 



The Weiss temperature 8 — 139 K and the broad max- 
imum around T max = 85 K, corresponding to ~61% of 
9, indicate sizable AF spin correlations. The magnetic 
ordering temperature TV = 28 K can be clearly traced 
by the divergence of the x(T) curves measured in low 
and high fields. While the low-field measurements show 
only an inconspicuous bend around T/y, the measure- 
ments above 2 T reveal a well-defined cusp. This diver- 
gence is due to the spin-flop transition in the AF-ordered 
phase. The spin-flop transition indeed takes place at 
about 1.7 T, as shown by the low-temperature magne- 
tization curve. 25 

Room-temperature ESR measurements yield a narrow 
line that can be fitted with a single powder-averaged 
Lorentzian. 25 Within resolution, the resonance field ap- 
pears isotropic and corresponds to the g value of 1.968 
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FIG. 3. (Color online) Top panel: Experimental magnetic 
susceptibility x{T) of Cr2[BPaOi2] measured at 5T (circles) 
and fit with the model of coupled bond-alternating S—3/2 
chains (line). Bottom left panel: Curie- Weiss fit to experi- 
mental \(T). Bottom right panel: Field dependence of x(T). 
Note the long-range AF ordering transition at Tn = 28K. 



which indicates weak spin-orbit coupling, as typical for 
Cr 3 +. 

The low-dimensional magnetism of Cr 2 [BP 3 0i2] is 
consistent with the presence of C^Og dimers in the 
crystal structure. However, isolated spin dimers should 
have a singlet GS at low temperatures. The fact that 
Cr 2 [BP 3 0i2] develops long-range AF order with a sizable 
Neel temperature of 28 K implies substantial interdimer 
couplings via the PO4 groups. While such couplings are 
abundant in vanadium and copper phosphates, 26-28 their 
identification is by no means a simple task. For exam- 
ple, in Cr 2 [BP 3 0i2] the shortest interdimer Cr-Cr dis- 
tances are 4.55 A along [0001] and 4.67 A in the (0001) 
plane, and it is impossible to decide a priori whether 
one of these pathways is more efficient, or both should 
be treated on equal footing. To develop a reliable mag- 
netic model of Cr2[BP 3 0i2], we perform extensive DFT 
calculations followed by QMC simulations, and provide 
detailed microscopic insight into the magnetism of this 
compound. 



IV. DFT CALCULATIONS 

To unravel the microscopic origin of the observed low- 
dimensional magnetic behavior, we apply DFT calcu- 
lations and evaluate the individual exchange couplings. 
Nonmagnetic (spin-unpolarized) calculations within both 
LDA and GGA yield a well-structured valence band 
(Fig. 4, top) with the total width of ~10 eV, crossed 
by the Fermi level ep at zero energy. This metal- 
lic electronic GS, contrasting with the green color of 
Cr2[BP 3 0i2], originates from the underestimation of the 



strong electron-electron repulsion within the Cr 3d shell. 
Technically, the band gap can be readily restored in a 
spin-polarized calculation, due to the sizable exchange 
splitting typical for Cr 3+ . However, its value is too 
small (1.29 eV in LSDA, 1.78 eV in GGA) 25 to account 
for the green color (~2. 2-2.5 eV) of Cr 2 [BP 3 12 ]. Ac- 
counting for electronic correlations is challenging, since 
many-body effects can not be properly described within 
the one-electron approach of conventional DFT function- 
al. 

Typically, multiorbital correlated insulators are de- 
scribed by an extended Hubbard model that comprises 
the kinetic terms (electron transfer), the on-site and in- 
tersitc correlations (Coulomb repulsion), and the on-site 
exchange (Hund's exchange). However, the large Hilbert 
space of such models impedes even an approximate nu- 
merical solution that would establish a simple relation 
between the parameters of the Hubbard model and the 
resulting magnetic couplings. Therefore, we restrict our- 
selves to a qualitative analysis of electron transfers. 

For the quantitative evaluation of magnetic couplings, 
we account for correlation effects using the mean-field- 
like DFT+J7 approach or hybrid (DFT + exact ex- 
change) functionals. The comparison to the experiment 
demonstrates the good agreement between different com- 
putational approaches, and underscores the validity of 
the qualitative analysis based on the electron transfers. 
Thus, this qualitative model approach could also be used 
as a starting point for more involved Cr 3+ systems with 
nontrivial coupling pathways. 

A. Qualitative model approach 

Conventional DFT (LDA and GGA) functionals are 
known for their accurate description of electron transfer 
processes, thus kinetic terms of the model Hamiltonian 
can be evaluated directly from the LDA/GGA band dis- 
persions for the Cr 3d states. Typical for an octahedral 
environment, these states split into two manifolds: half- 
filled tig orbitals centered at Ef, and empty e g orbitals 
that lie ~1.7eV higher in energy owing to the sizable 
crystal-field splitting. 

Altogether, the ti g and e g manifolds comprise 20 
bands, in line with five 3d orbitals per Cr atom and four 
Cr atoms in a unit cell. To evaluate the transfer inte- 
grals t™- m (m and m! are orbital indices, % and j are site 
indices), we map these 20 bands onto an effective five- 
orbital tight-binding model, and parameterize this model 
using the Wannier functions (WF) technique. 29 In the 
WF basis (Fig. 5), the couplings f^ m are evaluated as 
nondiagonal matrix elements. The resulting t™- m show 
excellent agreement between the tight-binding model and 
the computed LDA/GGA band dispersions (Fig. 4, bot- 
tom). 

Sizable transfer integrals are found for four Cr-Cr 
pathways (Table III): the intradimer transfer (ti), the 
interdimer tranfer along c (i^), the shortest interdimer 
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r a 
k vector 

FIG. 4. (Color online) Electronic structure of Cr2[BPsOi2]. 
Top: LDA density of states. The Fermi level £f is at zero 
energy. The contribution from Cr states is shown by shaded 
filling. Bottom: LDA bands in the vicinity of the Fermi level 
and the fit with the tight-binding model based on Wannier 
functions. 



pathway in the ab plane (£^2), and the longer interdimer 
pathway along the [223] direction (ij c i)> see Fig. 1. Here, 
il x and ti C 2 run via triple bridges of the PO4 polyhedra 
(along c), whereas t icl runs via a single PO4 tetrahedron. 
Other couplings are negligible, as they involve the trans- 
fer through at least two PO4 tetrahedra in a row. 

Although an explicit expression relating exchange in- 
tegrals Jij to transfer integrals t 7 ^" 1 is presently not 

available, the qualitative comparison of t™ m terms elu- 
cidates different contributions to the magnetic exchange. 
According to the Goodenough-Kanamori rules, the hop- 
pings between the half-filled states (both m and m' be- 
long to the tig subspace) are responsible for AF cou- 
plings, whereas the hoppings between the half-filled and 
empty states (m belongs to t 2g , m! belong to e g ) give rise 
to ferromagnetic (FM) interactions. 

First, we consider the intradimer t™ m hoppings, where 
the couplings between different t 2g orbitals are dominant: 
^xy.xz _j.xz,yz = tv z ' xy = — 125meV. At the same time, the 
largest hopping between the ti g and e g orbitals is smaller 

(t^ y ' z ~ r = — 97meV), thus hinting at the AF nature 
of Ji. For the interdimer coupling t' x the difference is 
even more pronounced: the hoppings between t2 g orbitals 
amount to — 124 meV, thus dominating over the hopping 
to the empty z 2 — r 2 orbital. Therefore, the interdimer 
exchange J[ should be also AF. 

The interchain coupling J- lC 2 is realized primarily via 



TABLE III. Leading transfer integrals t™ m (notation ac- 
cording to Fig. 1), where m and m' are orbital indices from 
the following set: \xy), \xz), \yz), \z 2 — r 2 ), and \x 2 ~y 2 ). All 
values are given in meV. For clarity, only one of two symmet- 
rically equivalent terms t™ m and m is shown. The Cr-Cr 
distances (do-Cr) are given for the 4K structure (Table I). 
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hoppings, hence the AF contribution is 
again dominant. In contrast, for the interchain coupling 
Jici, the t2 g — > t2g hoppings are close to zero, whereas 
the t2 g — > e g hoppings are still sizable. Therefore, the 
overall exchange should be FM. 

The above qualitative analysis is confirmed by elabo- 
rate DFT+C7 and hybrid-functional calculations reported 
below. More importantly, the proposed couplings are in 
agreement with the experimental magnetic structure that 
features antiparallel spins along J\, J[, and J; C 2 (Fig. 1, 
bottom left panel) . The respective AFM couplings estab- 
lish parallel spins along the J lc i bonds that are indeed 
FM. Therefore, the spin lattice of Cr2[BP30i2] is non- 



\xy) \xz) \yz) \z 2 — r 2 ) \x 2 — y 2 ) 

FIG. 5. (Color online) Cr2[BP30i2]: Wannier functions for five Cr 3d orbitals (see subscripts). Cr, P and O atoms are depicted 
by large (red), intermediate-size (green) and small (red) spheres, respectively. For each of the five Wannier functions, a Cr2 0g 
dimer and three PO4 tetrahedra are shown. The projection is similar to the bottom right panel of Fig. 1. 



frustrated. Because J{ 3> Jici7^ic2, this spin lattice can 
be considered as quasi- ID, with bond-alternating Ji — J[ 
chains running along the c direction. 



B. Total-energy calculations 

A sizable exchange splitting typical for Cr 3+ com- 
pounds readily opens a band gap (albeit underestimating 
the experimental value, see Ref. 25) in the spin-polarized 
calculations, so that the magnetic energy can be eval- 
uated directly from LSDA or GGA total-energy calcu- 
lations. In particular, total energies E tot correspond- 
ing to different collinear magnetic configurations can be 
mapped onto classical Heisenberg model, thus yielding 
the exchange integrals J if. 

Etot = Eq + E magn = Eq + JijSi ■ Sj. (3) 

(id) 

We apply this approach to evaluate the intrachain and 
interchain couplings. In LSDA, we obtain J\ = 151 K 
and J[ = 82 K as well as the FM J icl = -8 K and the AF 
Jic2 = 2 K. Spin-polarized GGA yields marginally smaller- 
values (Table IV). To challenge these estimates, we ad- 
dress the Weiss temperature which is a linear combina- 
tion of magnetic couplings. 30 For the Cr 2 [BP 3 0i2] spin 
lattice, the expected Weiss temperature is: 

6 = \S{S + l)(.h + J[+ 6Jici + 3J ic2 ), (4) 

where each coupling is summed up according to its co- 
ordination number in the spin lattice. Adopting the 
values of J* from LSDA (GGA), we obtain 6 -242K 
(202 K), which is 50-70% larger than the experimental 
value 6> = 139K. 

This discrepancy primarily originates from the poor 
description of electronic correlations within LDA or 
GGA, and calls for the application of more elaborate 



computational approaches. For strongly correlated in- 
sulators, such as Cr2[BP30i2], the mean-field DFT+C/ 
method is a natural choice. An intrinsic problem of this 
method is the double counting (DC) of the correlation 
energy already present in LDA or GGA. This correla- 
tion energy should be subtracted from the total energy 
of the system. In the widely used DC corrections, de- 
noted around-mean-field (AMF) and fully-localized limit 
(FLL), 31 the subtracted energy corresponds to the energy 
of the uniformly occupied state or the state with inte- 
ger occupation numbers, respectively. Previous DFT+C/ 
studies on Cr 3+ materials 32 did not render any of the 
two schemes preferable, thus we apply both and compare 
their results with the experiment. 

Besides the DC, the DFT+C/ method introduces two 
free parameters: the on-site repulsion Ud and the on-site 
exchange Jd- While an empirical evidence favors Jd — 
1 eV for 3d elements, the values of Ud can substantially 
vary depending on the electronic configuration and local 
environment of the magnetic ion. For Cr 3+ , we varied 
Ud in the range between 2 and 4eV. 32 ' 33 

The resulting DFT+C/-based exchange integrals are 
listed in Table IV. We find that the exchange integrals 
are weakly dependent on the DFT exchange-correlation 
potential (LSDA or GGA). In contrast, the DC correc- 
tion plays a more substantial role: the calculations within 
AMF yield considerably larger couplings than the calcu- 
lations within FLL for the same value of Ud- Besides, J; c i 
turns out to be AF in AMF, but FM in FLL. This con- 
spicuous difference between different flavors of DFT+C/ 
necessitates an additional examination of this problem 
by an independent technique. 

To this end, we resort to the HSE06 hybrid functional 
that is free from the DC problem. Here, the DFT ex- 
change is mixed with a fraction of exact (Hartree-Fock) 
exchange. 34 The HSE06-based exchange integrals (Ta- 
ble IV) are very similar to those obtained using DFT+C/ 
FLL with U d = 3-0 eV. To check the stability of the HSE06 
results, we varied the parameter /3, which reflects the ad- 
mixture of the exact exchange to the standard DFT ex- 



change (/3 = 0.25 in the original definition of HSE06). As 
follows from the resulting values (Table IV), the increase 
in j3 from 0.2 to 0.3 is accompanied by the ~25% reduc- 
tion in the exchange couplings, while the J[/Ji ratio is 
essentially unchanged. Therefore, the effects of increas- 
ing (3 in HSE06 and increasing Ud in DFT+U are some- 
what similar, despite the disparate physical meaning of 
the Ud and /3 parameters. 

On a qualitative level, DFT+U FLL, hybrid func- 
tional, and model approach concur with each other on 
the nature of magnetic couplings in Cr2[BP30i2]- By 
contrast, DFT+U AMF predicts the different sign of 
Jici, which would render the spin lattice weakly frus- 
trated. However, this scenario looks very unlikely be- 
cause the independent model analysis (Sec. IV A) yields 
FM J icl . Additionally, DFT+U AMF systematically 
overestimates the absolute values of Ji and J[ (see ex- 
perimental values in Sec. V) and should probably be dis- 
couraged in the case of Cr 2 [BP 3 0i2]. 35 

Apart from the overall energy scale of J\ and J[ that 
is dependent on Ud or j3 (Table IV), different flavors of 
DFT+C/ as well as HSE06 predict different J' 1 /J 1 ratios. 
Our computational results yield J[/Ji = 0.39 — 0.94, so 
that J[ is smaller than Ji but certainly large enough to 
ensure the quasi- ID nature of the spin lattice. Unfortu- 
nately, the more precise evaluation of the J[/Ji ratio lies 
beyond the capabilities of present-day DFT-based tech- 
niques and should be addressed by numerical simulations 
of the spin model and subsequent fitting to the experi- 
mental data (see Sec. V). 

Finally, we take into account the spin-orbit (SO) cou- 
pling, and quantify the single-ion anisotropy D using the 
DFT+U+SO calculations within the FLL using vasp. 
To this end, we calculate total energies (Fig. 1, bottom 
left) of the AF GS with the magnetic moments aligned 
parallel and perpendicular to the hexagonal axis (c-axis). 
For Ud = 3 eV, the energy difference between the two con- 
figurations is only 1.3 K, which is in accord with the low 
spin-flop field of- 1.7 T. 

In summary, our band structure calculations arrive at 
a scenario of Heisenberg spin chains with two alternat- 
ing nearest-neighbor AF interactions J\ and J[. The 
magnetic chains are weakly coupled in a non-frustrated 
manner by the FM J lc i and the AF J ic2 . While both 
interchain couplings are apparently weak (below 2 K in 
terms of the absolute value), the energy scale of the intra- 
chain couplings, as well as the precise alternation ratio, 
are rather sensitive to details of the computational proce- 
dure and require further refinement applying simulation 
techniques to fit the experimental data. 



V. QMC SIMULATIONS AND COMPARISON 
WITH THE EXPERIMENTS 

To refine the parameters of the microscopic magnetic 
model, we simulate the temperature dependence of the 
reduced magnetic susceptibility x* f° r various ratios of 



TABLE IV. Leading exchange integrals Ji (in K) calculated 
with different DFT-based methods. For the DFT+U results, 
the double-counting correction (DCC) and Ud value (in eV) 
are given. For the HSE06 hybrid functional, the admixture of 
the exact exchange (/3) is specified. 



£ xc (DCC) 


U d or/3 


Ji 


Ji 


Jic2 


Jicl 


Ji/Ji 


LSDA 




150 


82 


2.5 


-7.8 


0.55 


GGA 




125 


63 


2.2 


-5.6 


0.51 




2.0 


109 


60 


2.3 


0.1 


0.62 


LSDA+U (AMF) 


3.0 


89 


55 


2.2 


1.7 


0.72 




4.0 


70 


50 


2.1 


2.8 


0.86 




2.0 


84 


49 


1.6 


-2.1 


0.58 


LSDA+U (FLL) 


3.0 


54 


38 


1.2 


-1.8 


0.71 




4.0 


33 


31 


0.9 


-1.6 


0.94 




2.0 


97 


48 


2.2 


1.0 


0.50 


GGA+U (AMF) 


3.0 


81 


44 


2.1 


2.3 


0.54 




4.0 


67 


41 


2.1 


3.3 


0.61 




2.0 


78 


40 


1.6 


-1.0 


0.51 


GGA+U (FLL) 


3.0 


54 


32 


1.2 


-1.0 


0.59 




4.0 


37 


26 


1.0 


-0.9 


0.69 




0.20 


66 


27 


1.2 


-2.5 


0.41 


HSE06 


0.25 


58 


23 


1.0 


-2.0 


0.40 




0.30 


51 


20 


0.9 


-1.7 


0.39 



J[/Ji, while keeping the interchain couplings constant 
Ji C 2 = — Jici=0.027i (see Table IV), and subsequently fit 
the simulated QMC curves to the experiment using the 
expression: 

Here, the fitted parameters are Ji, the Lande factor g, 
and the temperature-independent contribution xo- 36 The 
results are presented in Table V (columns 2 and 3). A 
general trend is the decrease in 3\ upon the increase in 
J[/Ji, so that the sum of the leading couplings J\ + J[ 
remains nearly constant (74.5 ± 2.0 K), see also Eq. (4). 
Still, smaller values of J[/Ji yield better agreement with 
the experimental curve, especially around the maximum 
in x{T). 

In general, additional restrictions for the model param- 
eters can be set by the fitted values of g. However, g is 
marginally dependent on J[/Ji and shows good agree- 
ment with the ESR value of 1.968 in the whole range 
studied. Therefore, simulation of further measurable 
quantities is needed to refine the model parameters. 

First, we address the magnetic ordering temperature. 
To this end, we calculate the spin stiffness p as a func- 
tion of temperature for finite lattices with up to N — 
L x x L y x L z = 13 824 spins. At the magnetic ordering 
transition temperature, the products piLi(i = x,y, z) are 
independent of Lj . Thus, Tn is the temperature at which 
the pL(T) curves of the different finite lattices cross. The 
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TABLE V. Results of QMC simulations and fitting to the 
experimental magnetic susceptibility. Ji and g are obtained 
from the fits to x(T). Tn/Ji is obtained from simulations of 
spin stiffness (see text). Sao is evaluated using Eq. (6). Tn 
and m are scaled using Ji and g from columns 2 and 3. 



J'l/Jl 


Ji (K) 


9 


Ts/Jx 


T N (K) 




m(fi B ) 


0.4 


51.84 


2.006 


0.584 


30.27 


1.128 


2.22 


0.5 


49.51 


2.002 


0.642 


31.79 


1.151 


2.27 


0.6 


47.10 


1.996 


0.689 


32.45 


1.164 


2.29 


0.7 


44.78 


1.993 


0.726 


32.51 


1.169 


2.30 


0.8 


42.58 


1.992 


0.755 


32.15 


1.165 


2.29 



difference between the simulated values of Tn and its ex- 
perimental value of 28 K is below 10% for the whole range 
of J'x/J\ (sec columns 4 and 5 of Table V). Similar to the 
fits of x(T), lower values of J[/Ji yield slightly better 
agreement with the experimental Tn = 28K. 

Another measured quantity, which can be used for 
comparisons between theory and experiment, is the or- 
dered magnetic moment to. In the classical 5 = 3/2 
Heisenberg model, the ordered magnetic moment to = 
gfJ-B>5 amounts to 3/ib, but the experimentally observed 
moment in Cr2[BP30i2] is substantially smaller (2.5 /xb, 
see Sec. MB). We estimate the reduction in to due to 
quantum fluctuations. To this end, we simulate the mag- 
netic structure factors S for finite lattices containing up 
to 2048 spins. 

For the propagation vector k of the magnetically or- 
dered GS, the ordered magnetic moment m is evaluated 
using the finite-size scaling: 



3S(k) 
N 



= St 



Vn n 

m = g^BSoo, 



(6a) 
(6b) 



where the fitting parameters are 5oo, o\ and 02, while 
g= 1.968 is adopted from the ESR. 37 In this way, we ar- 
rive at the values that underestimate the experimental 
result by ~ 0.2 /xb (last column of Table V). Since this 
offset is nearly independent of J(/Ji, the interchain cou- 
plings might be the origin of this discrepancy. Indeed, 
increasing the interchain couplings by a factor of two 
(Jic2 = — -Act = 0.04 Ji) yields to = 2.49/x b , in excellent 
agreement with the experiment. However, larger inter- 
chain couplings substantially increase the Neel temper- 
ature and worsen the agreement between the simulated 
and experimental x(T). Since the experimental value of 
Tn is more precise than to, we argue that the parameters 
Ji C 2 = — J id =0.02 Ji are preferable, while the 0.2 /jb de- 
viation between the experimental and simulated values 
of to is still reasonable. Indeed, the experimental esti- 
mate of m relies upon the available magnetic form-factors 
as well as on the accuracy of the subtraction procedure. 
Therefore, we can not exclude a systematic experimental 
error in m that should be verified by neutron experiments 
on single crystals. 



As follows from Table V, different values of J[/Ji 
yield rather similar Tn and mo- This surprisingly ro- 
bust behavior originates from peculiar properties of spin- 
3/2 alternating Heisenberg chains that show low spin 
gaps in the 0.41 < J[/Ji < 1 range between the crit- 
ical points at J[/Ji = 0.41 and J[/J\ = 1, where the 
spin gap vanishes (see Sec. VI and Ref. 38). To eval- 
uate the optimal J'i/Ji, we trace the general trends. 
First, smaller J[/Ji ratios yield better agreement with 
the experimental Tsr. Second, too is nearly indepen- 
dent of J'xl J\ for J' x l J\ >0.5. Thus, the parametrization 
Ji : J[ : Jid ■ Jic2^1 : 0.5 : -0.02 : 0.02 with Ji~50K 
is the optimal choice, since it conforms to the experimen- 
tal x(^) dependence, accurately reproduces the experi- 
mental Tn as well as the g-factor, and is in reasonable 
agreement with the experimental too- 



VI. QUANTUM-CLASSICAL CROSSOVER 

For a correct application of classical models, it is 
crucial to know at which point quantum effects be- 
come relevant and the classical approximation breaks 
down. In this respect, the quasi-lD magnetic model of 
Cr 2 [BP 3 0i2] is a promising candidate. The leading ex- 
change couplings Ji and J[ form alternating chains, while 
the interchain couplings are substantially smaller. Be- 
sides, Cr2[BP30i2] lacks frustration, thus quantum MC 
simulations can be performed. 

First, we study nearest-neighbor spin correlations us- 
ing QMC and classical MC simulations. The difference 
between the two results is directly related to quantum 
fluctuations: while QMC correctly accounts for the quan- 
tum behavior at finite temperatures, the classical MC 
method captures thermal fluctuations, only. 

In the following, we consider a ID J\ — J[ magnetic 
model, where the interchain couplings are neglected com- 
pletely. For J[/ Ji = 0.5, we simulate the temperature de- 
pendence of the diagonal spin correlations (S*S* +1 ) and 
(5f +1 Sf +2 ), where i = corresponds to the stronger cou- 
pling J\ , while i = 1 corresponds to J[ . 

The resulting curves are plotted in Fig. 6. In the 
high-temperature range (see inset in Fig. 6), the classical 
and quantum results are practically indistinguishable. At 
T ~ Ji , the curves start to deviate significantly. The clas- 
sical MC curves exhibit an asymptotic behavior and join 
at T = 0. This is in line with a complete elimination of 
thermal fluctuations, thereby the spin correlations reach 
their extremal value for S — 3/2 spins, which amounts to 
-[5(5+ l)]/3 = -1.25. 

Although the QMC curves exhibit qualitatively simi- 
lar behavior at low temperatures, they saturate at sub- 
stantially higher values of (S*Sf +1 ). Dissimilar to the 
classical model, the i = ( J\) and i = 1 (J{) QMC curves 
saturate at different values of (SfSf +1 ), indicating the 
onset of dimerization. 

In the vicinity of the broad maximum in the magnetic 
susceptibility (T max ~2 Ji), and even at lower tempera- 
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FIG. 6. (Color online) Temperature dependence of the 
nearest-neighbor diagonal spin correlations (5'f5'f +1 } in the 
S — 3/2 alternating Heisenberg chain model (Ji/Ji = 0.5). 
QMC results are plotted with squares (Ji) and circles (J[), 
the classical MC results are depicted with pluses (Ji) and 
crosses (J[). The dashed lines denote the ED results for the 
quantum model. The ^S(S + 1) line marks the maximal di- 
agonal correlation for S= 3/2 spins. Inset: at high temper- 
atures (T>2Ji), the QMC and MC curves become almost 
indistinguishable. 



tuxes, the classical MC simulations are in excellent agree- 
ment with the quantum model. Therefore, experimental 
thermodynamic data [e.g., x(^) an(1 magnetic specific 
heat] for T>J\ should be well reproduced within the 
classical approximation. This empirical rule follows ear- 
lier experimental and numerical results for uniform spin- 
3/2 chains. 9 It can also be generalized to quasi-2D and 3D 
systems, owing to further suppression of quantum fluc- 
tuations in higher dimensions. However, in 0D systems, 
such as isolated dimers, quantum effects are expected to 
play a more substantial role. 

To explore the evolution of quantum fluctuations by 
a gradual crossover from ID to 0D regime, we evalu- 
ate the GS diagonal spin correlations for different ratios 
J'x/Jx- Although QMC simulations are not applicable at 
T = 0, the temperature evolution of spin correlations in 
Fig. 6 suggests a marginal change, only, of (SfS? +1 ) be- 
low 0.03 Ji. Thus, finite-temperature spin correlations 
simulated at T = 0.03 Ji can be regarded as a rather ac- 
curate measure of GS (T = 0) spin correlations. 

Fig. 7 (top) shows the resulting {SfSf +1 )(J[/Ji) de- 
pendencies. In the uniform-chain limit (J[/J\ = 1), 
the correlations for i = and 1 coincide. Dimerization 
gives rise to a quasilinear behavior observed down to 
J'x/Jx — 0.4. At this point, both curves exhibit a pro- 
nounced kink, and their further evolution is different: 
the i — (Ji) curve rapidly drops almost down to the 
classical value \S[S + 1], while the i = 1 (J[) curve ex- 
hibits a quasilinear growth up to zero correlation in the 
J'xl J \ = limit. This behavior is remarkably different 
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FIG. 7. (Color online) Top: nearest-neighbor diagonal spin 
correlations {SfSf + i) in the S = 3/2 alternating Heisenberg 
chain model as a function of the alternation ratio J{/Ji at 
T = 0.03Ji. Dashed lines are guides to the eye. Bottom: root 
mean square of Oi = ({S*S* +1 ) — + 1]) (i = 1, 2), reflect- 

ing the discrepancy between the quantum and the classical 
models. 



from the classical case, where the constant spin-spin cor- 
relation of — + 1) is expected. 

To sharpen the crossover at J'xl Ji~0.4, we 
consider root mean square deviation for ui = 
((S?Sf +1 ) - lS{S + 1)) (i = l, 2), which reflects the dif- 
ference between the quantum and classical models. The 
resulting dependence is shown in Fig. 7 (bottom). For 
0.5 < J'\l J \ < 1, only a marginal increase is observed, 
while at J'xl J \ — 0.4 the slope changes distinctly. For 
J'xl Jx <0.4, a quasilinear behavior is restored. 

The drastic change in the spin-spin correlations around 
J'x/Ji — 0.4 indicates the critical point of bond- 
alternating spin-3/2 chains. While such chains generally 
have gapped GS, they become gapless in the uniform- 
chain limit {J[/Jx — 1) and at the critical point of 
J'xl J \ — 0.41, according to density-matrix renormaliza- 
tion group 38 and QMC studies. 39 Qualitatively, the criti- 
cal point at J'xl J \ — 0.41 can be understood as the transi- 
tion between two different valence-bond-solid (VBS)-type 
GSs. 

The VBS state entails singlet pairs formed by indi- 
vidual spin- 1/2 entities comprising spin-3/2 sites of the 
lattice. 40 Although it is not an exact GS of the Heisen- 
berg Hamiltonian in Eq. (1), the references to relevant 
VBS states are instrumental in understanding peculiari- 
ties of quantum magnets. The spin-3/2 chain can form 
two VBS states, the (3,0) state having three singlet pairs 
on one bond and no pairs on the contiguous bonds, or 
the (2,1) state having two singlet pairs on one bond and 
one singlet pair on each of the neighboring bonds. At low 
J'xl 1 o h, alternating spin-3/2 chains are close to the (3,0) 
regime, as shown by the drastic difference between the 
spin-spin correlations for Ji and J[. At J'xl J \ > 0.41, 
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the system resembles the (2,1) VBS 39 and eventually de- 
velops the GS of a uniform spin-3/2 chain with equal cor- 
relations on the Ji and J[ bonds (Fig. 7). Note that the 
low-temperature physics of alternating spin-3/2 chains is 
essentially governed by quantum effects and can not be 
reproduced in the classical approximation. 



VII. DISCUSSION AND SUMMARY 

The formation of alternating spin chains contrasts with 
the presence of isolated Cr 2 Og dimers in the crystal struc- 
ture of Cr2[BP30i2]. However, such differences between 
the structure and the ensuing spin lattice are not un- 
usual, especially for phosphates, vanadates, and other 
compounds with complex anions that are capable of me- 
diating the superexchange. For example, a large su- 
perexchange coupling via the PO4 tetrahedra between 
the structural dimers is a salient feature of the spin- 
1/2 alternating-chain magnet (VO) 2 [P 2 7 ] (Ref. 26). In 
spin-1/2 V 4+ and Cu 2+ compounds, intcrdimcr couplings 
often exceed the intradimer exchange, whereas the latter 
is small or even negligible, as in CU2P2O7] (Ref. 28) or 
VO[HP0 4 ] ■ §H 2 (Ref. 41). 

The Ci - 2[BP30i2] case is different, though, because the 
intradimer coupling J\ exceeds the interdimer coupling 
J[ . The reason behind the sizable intradimer coupling is 
the strong direct exchange that was previously observed 
in LiCr0 2 . 32 While the Cr-O-Cr angles amount to 87.8° 
and impede the AF superexchange, the direct overlap of 
Cr d orbitals is strong enough to facilitate the sizable 
AF Ji. This mechanism is not operative in typical spin- 
1/2 systems with structural dimers, because V 4+ cations 
form VO5 square pyramids and remain on opposite sides 
of the basal plane of such pyramids, thus featuring only 
a weak direct overlap of the magnetic d xy orbitals. In 
the case of Cu 2+ , the magnetic orbital has the d x 2_ y 2 
symmetry, unfavorable for the direct overlap. 

It is also instructive to compare different interdimer 
couplings in Cr2[BP30i2]. Similar to other transition- 
metal phosphates, 27,28 the Cr-based Wannier functions 
show sizable contributions of Cr 3d and O 2p states, only 
(Fig. 5). Phosphorous states weakly contribute to the 
magnetic orbitals and play a minor role in the superex- 
change running via the Cr-0...0-Cr pathways. The ef- 
ficiency of these pathways intimately depends on details 
of the crystal structure. 

Remarkably, Cr2[BP30i2] not only entails the spin lat- 
tice of alternating chains, but also features the alterna- 
tion ratio of J[/Ji — 0.5, which is close to the criti- 
cal point of the gapless ground state at J[/Ji — 0.41. 
The proximity to this critical point is one of the rea- 
sons behind the long-range AF ordering and the siz- 
able Neel temperature in Ci'2[BP 3 0i2]. Spin-1/2 sys- 
tems with the larger J[/J\ — 0.6 (i.e., weaker dimer- 
ization) and comparable interchain couplings 42 still have 
the spin-gap ground state, as in Pb 2 [V 3 09] (Ref. 43) 
and Ag(VO)[As0 4 ] (Ref. 44). The large number of in- 



terchain couplings per magnetic site (the coordination 
number amounts to six and three for J; c i and Jj c2 , re- 
spectively), as found in Cr 2 [BP 3 0i 2 ], may also reduce 
quantum fluctuations, thus facilitating the long-range AF 
order in Cr2[BP30i2]. By contrast, an isolated alternat- 
ing S = 3/2 chain with J[/Ji — 0.5 would feature an 
excitation gap of ^0.25 J\ with no long-range order. 39 

While the critical point at J[/Ji — 0.41 has been 
widely studied theoretically, 38 ' 39 experimental probes of 
this regime and even experimental examples of spin-3/2 
alternating chains have not been reported. Therefore, 
Ci"2[BP30i2] may be interesting as an experimentally 
available spin-3/2 alternating-chain system in the vicin- 
ity of this critical point. Although Cr2[BP30i 2 ] is not 
perfectly ID and features the long-range AF order owing 
to non-zero interchain couplings, an experimental study 
of spin-spin correlations with, e.g., inelastic neutron scat- 
tering, could be instructive. Moreover, the application 
of external pressure might change the J[/ Ji ratio, thus 
giving access to the peculiar evolution of spin-spin corre- 
lations across the critical point (Fig. 7). Detailed studies 
of the isostructural phase Fe2[BP30i2] (Ref. 23) may be 
instructive as well, because this compound should fea- 
ture a similar spin lattice in the more classical regime of 
spin- 5/2. 

In summary, we extensively characterized the magnetic 
behavior of Cr2[BP30i2] using magnetic susceptibility, 
neutron diffraction, and ESR measurements. The long- 
range AF order established below Tn = 28 K is described 
with the propagation vector k = 0, whereby the spins 
on nearest-neighbor Cr atoms within the ab plane as well 
as along the c-axis direction are antiparallel. On the 
microscopic level, Cr2[BP30i2] features S — 3/2 Heisen- 
berg chains with an alternation of the nearest neighbor 
couplings Ji and J[. The ratio J[/Ji — 0.5 renders 
Ci"2[BP30i2] as an interesting model system lying close 
to the critical point at J[/Ji — 0.41, where the spin- 
3/2 alternating chain has the gapless ground state. The 
chains are coupled by two nonequivalent interchain ex- 
changes: the FM Jid and the AF J; C 2, both of the order 
of 1-2 K. The microscopic model is in excellent agreement 
with the experimental magnetic structure (see Fig. 1): 
the nearest-neighbor spins within a chain are antiparal- 
lel, in accord with the AF nature of J\ and J[. More- 
over, even the interchain coupling regime is compatible 
with FM J ic i and AF J ic2 , thus Cr 2 [BP 3 0i2] lacks any 
appreciable frustration effects. 
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FIG. 8. Field dependence of the magnetization of CT2P3P3O12] measured at 2 K. The lines 
are guide-for-the-eye. The inset shows the held derivative of the magnetization. 
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FIG. 9. Electron spin resonance spectrum (circles; first derivative of absorbed microwave 
power dP/dB) taken at 9.5 GHz (X-band). Solid line denotes a powder-averaged Lorentzian 
line with g ± » ^ = 1.9680 ± 0.0005 and linewidth AB ± = 8.6 ± 0.2 mT and ABy = 
18.8 ± 0.2 mT. 
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FIG. 10. Cr 2 [BP 3 0i 2 ]: GGA (solid line) and LSDA (dashed line) density of states for the 
ferromagnetic solution. Positive and negative values refer to the majority and minority spin 
channel, respectively. Note the different values of the band gap: — 1.78 eV in GGA and 
-1.29 eV in LSDA. 



